[Analog_Engine.git] / scripts / TP1 10.29 Two masses fixed to walls.LACE

# TP1 10.29 Two Masses fixed to walls
# Two equal masses m are connected to each other by a spring (spring constant D) and on either side to a fixed wall.

# m*s1'' + 2D*s1 - D*s2 = 0 => s1'' = -D/m*s1 -D/m*s1 + D/m*s2
# m*s2'' + 2D*s1 - D*s1 = 0 => s2'' = -D/m*s2 -D/m*s2 + D/m*s1

coefficient(1): D/m_1 # D/m for s1
coefficient(2): D/m_2 # D/m for s2, identical to D/m_1
coefficient(3): -1 -> -s1_0 # initial position of s1
coefficient(4): -1 -> -s2_0 # initial position of s2

iintegrate -D/m*s1, -D/m*s1, D/m*s2 -> -s1' # input is s1''
iintegrate -s1' -> s1
   IC: -s1_0
s1 * D/m_1 -> D/m*s1
invert D/m*s1 -> -D/m*s1

iintegrate -D/m*s2, -D/m*s2, D/m*s1 -> -s2' # input is s2''
iintegrate -s2' -> s2
   IC: -s2_0
s2 * D/m_2 -> D/m*s2
invert D/m*s2 -> -D/m*s2

output(x): s1
output(y): s2
Commit	Line	Data
56968557 P	1	# TP1 10.29 Two Masses fixed to walls
	2	# Two equal masses m are connected to each other by a spring (spring constant D) and on either side to a fixed wall.
	3
	4	# ms1'' + 2Ds1 - Ds2 = 0 => s1'' = -D/ms1 -D/ms1 + D/ms2
	5	# ms2'' + 2Ds1 - Ds1 = 0 => s2'' = -D/ms2 -D/ms2 + D/ms1
	6
	7	coefficient(1): D/m_1 # D/m for s1
	8	coefficient(2): D/m_2 # D/m for s2, identical to D/m_1
	9	coefficient(3): -1 -> -s1_0 # initial position of s1
	10	coefficient(4): -1 -> -s2_0 # initial position of s2
	11
	12	iintegrate -D/ms1, -D/ms1, D/m*s2 -> -s1' # input is s1''
	13	iintegrate -s1' -> s1
	14	IC: -s1_0
	15	s1 * D/m_1 -> D/m*s1
	16	invert D/ms1 -> -D/ms1
	17
	18	iintegrate -D/ms2, -D/ms2, D/m*s1 -> -s2' # input is s2''
	19	iintegrate -s2' -> s2
	20	IC: -s2_0
	21	s2 * D/m_2 -> D/m*s2
	22	invert D/ms2 -> -D/ms2
	23
	24	output(x): s1
	25	output(y): s2